oo 

o 

O 



Temperature dependent optical conductivity of undoped cuprates with weak exchange 
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The optical conductivity a-{u)) is calculated at finite temperature T for Cu02 chain clusters within 
a pd-Hubbard model. Data at T = 300 K for Li2Cu02 are reanalyzed within this approach. The 
relative weights of Zhang-Rice singlet and triplet charge excitations near 2.5 and 4 eV, respectively, 
depend strongly on T, and a rather dramatic dependence of a{Lo) on the ratio of the V* to 2"'' 
neighbor exchange integrals is predicted. On the basis of these results, information about exchange 
interactions for frustrated edge-shared cuprates can be obtained from T-dependent optical spectra. 
Our results are also relevant for magnetically weakly coupled wide-gap insulators in general. 

PACS numbers: 71.15.Mb, 71.27.-|-a, 74.72. Jt, 75.10. Jm, 75.10.Pq, 78.20.-e 



in 



I 

O 

o 



> 

en 
m 

00 

p 
o 

00 

O 



X 



Standard wisdom on wide gap insulators says that their 
optical spectra above 1 eV excitation energy are hardly 
affected by temperature T, magnetic fields H, and by 
the magnetic nature of their ground state (GS). More- 
over, spin and charge degrees of freedom are often de- 
coupled in one dimension (ID). Here, we present exact 
theoretical results qualitatively valid for several cuprates 
to be specified below which prove just the opposite in all 
these respects. Our study bases on the fact that different 
magnetic states with different symmetries obey different 
selection rules, which - in case of soft magnetic materi- 
als - can lead to sizable T-dependence of the optical spec- 
tra. For systems with small enough exchange integrals, 
the energy difference between the GS and the excited 
magnetic states of the system will be small enough so 
that they can be thermally populated. This will cause a 
strong T-dependence of various response functions, here 
presented for the case of the optical conductivity <y{Lu). 
For an illustration of our approach we consider Li2Cu02 
being of current interest^'^'^'^'^'^'^'^'^ i^°i^^'^^ , structurally 
simple, and most importantly, where single crystal data 
are available in a broad w-range^^. This system stands 
for a class of frustrated spin- 1/2 chain materials with 
both small ferromagnetic (FM) NN Cu-Cu exchange cou- 
pling Ji and AFM next-nearest neighbor (NNN) Cu-Cu 
in-chain exchange J^r-i^ in terms of a ID spin-1/2 Heisen- 
berg model 



Hs 






Jl^i^i+l + J2^i^i+2 + Js^i^i+S + 



(1) 



In Li2Cu02 chains running along the 6-axis are formed 
by the edge-sharing of Cu04 plaquettes. The most ex- 
citing puzzle addressed here is the missing of the Zhang- 
Rice singlet peak (ZRS) in reflectivity and electron en- 
ergy loss spectroscopy (EELS) data at T = 300 K^. Its 
detection in resonant inelastic x-ray scattering (RIXS) 
spectra is under debate, too^. There is also no con- 
sensus on the role of two scenarios for the FM in-chain 
order below the Neel temperature Tn = 9 K. Scenario I 
is given by a dominant FM Ji— defined by the inequal- 
ity a = —J2/J1 < ac w 0.25, if J3 is weaki^, where 



ac marks the transition from a high-spin FM to a low- 
spin spiral- like ground state (GS). While in scenario II 
it stems from a specific AFM inter-chain coupling and a 
single chain would be at a > aj^^^^. 

A first microscopic theoretical description^ of the mag- 
netic susceptibility xC^) ^md of a{uj) at zero temperature 
on the basis of a Cu Sd^y O 2px,y Hubbard model (a; and 
y along the b- and c-axes, respectively, see Figs. 1,2 of 
Ref. (ij) predicts charge transfer excitation of a Cu hole 
into a ZRS near Wzns= 2.25 eV. Its missing observation 
(at T— 300 K) was ascribed to poor resolution^, or to 
uncoupled CUO4 unitsi^ We show, that a ZRS visible 
at low T would be strongly suppressed at 300 K. How- 
ever, if the Hamiltonian parameters are chosen to im- 
prove the description of O Is x-ray absorption (XAS)i, 
optical, EELSii^, and RIXS data^, then the ZRS absent 
at r = 0, appears with rising T but is accompanied by a 
Zhang-Rice triplet (ZRT) contribution at lo^rt ~ 4 eV in 
(j{uj) (superimposed with further transitions). The (j(uj) 
data are well described this way. 

The used Hubbard Hamiltonian reads (cf. also Ref. [l|) 



Wh = 'Wk + Wc 



K 



(2) 



Its kinetic. Coulomb, and exchange contributions are 

i i.j,s 

He = ^ UiUi^riii + 2 X! ^»J"«"J' '^-^ ^^^nis, (4) 

i i^j s 



i^j,ss' 



where i and j run over all GvL-2)dxy and 0-2px^y orbitals 



and s is the spin index. Except for the Si the Hamilto- 
nian parameters are the same as in Ref. |l|: Ud =8.5 eV, 

^P:i,Pv = Up — 

2.9 eV for the O onsite inter-orbital repulsion, 
0.6 eV is the FM Hund's rule coupling. 



Up =4.1 eV for the intra-orbital and 

2Kp^ 

where K, 



Pi:Pv 



As in Rcf. 1 isotropy of the FM Cu-0 exchange inte- 
grals Kp^ = Kp^ci = Kp^d = 0.05 eV was used for the 
sake of simphcity. Polarized O Is XAS measurements^ 
with the electric field vector in x and y direction, respec- 
tively, indicated a nearly isotropic O 2p-hole distribution 
in the xy-plane. This indicates a condition Up^ ~ Up 
within XAS error bars for the O 2p hole occupation 
numbers in the GS. In order to achieve this despite the 
anisotropic CUO4 plaquette geometry, Sp^-Sp^ = 0.2 eV 
was taken as distinct from Ref. 1 where Sp^ = Sp was 
assumed. In order to reproduce the first strong peak in 
17(0;) near 4.4± 0.2 eV (Figs. 1,2), the mean O onsite 
energy Apd = {Sp^ + £pj,)/2 — Sd has additionally been 
up-shifted by 0.5 to 3.7 eV. To demonstrate the strikingly 
distinct a^uj) caused by these moderate changes, calcu- 
lations were performed with the £i-values of Ref. |l|, too. 
Hereafter these two models (2) are called Ml— and M2 
(above choice). We will show that they lead to different 
magnetic GS yielding different T-behavior of cr(uj). 

Naturally, the dominating exchange integrals Ji and J2 
determine the spectrum of low-energy excited states |i/) 
(spin excitations) of the spin model (1). We found the J- 
values from projecting the Hamiltonian (2) onto (1). At 
a given T and possibly in the presence of a magnetic field 
H, the optical conductivity (j{to, H, T) of (2) is obtained 
from the (t^{uj) with the initial spin states | v) of (2)^^: 

(j{uj,T) = ^w,{T)[l-cxp{~huj/k^T)]<j,{uj),{6) 



Y.U' 9i^' expi-E^'/ksT) ' 



(7) 



where for H = {) the spin degeneracy is gi, = 2Si, + 1 . For 
instance, for the largest cluster, CU6O14, which can still 
be handled by the exact diagonalization of the Hamilto- 
nian (2) there are 2^—64 low-energy spin states: 5 x (S* = 
0), 9 X (5 = 1), 5 X (S* == 2), and 1 x (S* = 3) multiplets. 
Due to the large optical transition energies Tiuj ^ k^T 
the thermal occupation of the final states in (6) can be 
ignored. In case of an applied external magnetic field H 
the S y^ states with finite Sz are Zeeman split which 
affects the Boltzmann probability Wi, = Wv (T, H) to find 
a cluster in a given spin state (5,^, Sz.v)- In this case the 
gi, are replaced by 



7^(i7,r) = 1 + 2^ cosh 



ly'gi^fisH 



k^T 



j.{H,T) = 2 2 cosh 



i2iy' - l)g^fi^H 



2k^T 



(8) 



(9) 



for even and odd chains, respectively, where (7l denotes 
the Lande-factor. All response functions reported be- 
low were calculated using exact diagonalizations and the 
common continued fraction method for Cu„02n-i-2 clus- 
ters. The (5-functions of the calculated a{uj)-spectra are 
convoluted with a Lorentzian broadening of 7^ ~ 0.35 eV 
at half width to compare them with optical data. 



Generally, the calculated a{uj) for Cu„02n-i-2 chain 
clusters at i/ = T = exhibit a multiple peak struc- 
ture: two well-pronounced peaks near 4 and 5 to 5.5 eV 
as shown in Fig. 1. Remarkably, there are marked dif- 
ferences in (t(w) between 2 and 4 eV depending whether 
the GS is low-spin (i.e. a spin spiral) or high spin (i.e. 
FM). For Li2Cu02, due to its closeness to the critical 
point ac (i.e. the FM-spiral transition), the issue de- 
pends sensitively on the J-values governed by the details 
of the parameters e^. A main issue of our analysis is that 
while the e^ values of Ml lead to a low-spin GS (spiral, 
a > ac) our choice M2 results in a high-spin GS (FM, 
a < Oc). (^i^) for a Cu60i4 chain with assumption of a 
high-spin GS (S=3) or a low-spin GS (S=0,1) are shown 
on Fig. 2. The height of the feature between 2 and 3 eV 
depends on the cluster size and on the total spin of the 
GS. Already a single CUO4 unit shows a double peak 
structure derived from transitions between an even par- 
ity singly occupied hybridized Cu 3d 02p state and odd 
parity non-bonding O states. Like the GS, all states are 
spin doublets. Since these transition energies do almost 
coincide with the mentioned above "high-energy" peaks 
for larger clusters, we regard them as intra-plaquette ex- 
citations. Inter-plaquette transitions may occur starting 
from dimers n = 2. If their GS is an 5 = state or a 
low-spin state for larger clusters (e.g. a doublet for odd 
chains), a peak appears between 2 and 3 eV. This tran- 
sition is forbidden, if the GS is the maximal spin state: 
i.e. an 5 = n/2 FM state. Since in the optical transi- 
tion S and Sz are conserved, the final state for an ex- 
cited dimer in the former case is again a singlet (low-spin 
state) with two holes at one of the two CUO4 units, one 
sitting mainly on Cu and one sitting mainly on O (simi- 
lar to a ZRS state on a CUO4 ^ unit). This transition is 
usually denoted as a ZRS excitation. In the FM case the 
excited state may contain a triplet (ZRT) of two holes 
on one plaquette, which in the optical excitation process 
occurs slightly below the main peak near 4 eV— . The 
intensities of the ZRS and ZRT transitions exhibit some 
positive finite size effect (x 1/n reflecting mainly the num- 
ber of available plaquettes for an inter-plaquette transi- 
tion. For a dimer with Ji < the GS is always a triplet 
since J2 is involved for n > 3, only. A detailed finite- 
size analysis will be given elsewhere. The ZRT-energy « 
4 eV found above is in excellent agreement with 4.1 eV 
reported in a recent RIXS study^ but in sharp contrast 
with the assignment given in Ref. |j where the RIXS peak 
near 2.1 eV has been ascribed to the ZRS using a low- 
spin GS. Probably, the 2.1 eV-peak should be ascribed 
to d-d excitationsiSiii. A similar feature has been found 
in CuGe03 (in (t{uj) at lower ui « 1.7 ± 0.1 eV) while 
the ZRS is clearly observed at 3.4 ± 0.2 eV-^^^. 

Exact diagonalization of (2) yields excitation energies 
tJzns = 2.7 eV, uJzRT = 4 eV for M2 and Wzrs = 2.25 eV, 
ujzRT = 3.7 eV for Ml. For both models this is consistent 
with Azn w 1.3 eV obtained as the distance between the 
lowest 5 = and the S = 1 levels in a single Cu04r^ 
plaquette with two holes. The a{u!,T) obtained from (7) 
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FIG. 1: (Color). Finite size eflFect for a{Lu). Left (right): 
GS given by the lowest (highest) spin state. The curves for 
various clusters are shifted by (n-l)xlOO(r2cm)~ . 



for both models are shown on Fig. 3, upper two pan- 
els. At 300 K Ml and M2 yield the same qualitative 
behavior: the ZRS is largely suppressed. The main peak 
position of the experiment at > 4eV is, however, much 
better reproduced by M2. We studied also the field de- 
pendence of a{uj,T, H) (Fig. 3, lower panel) and found 
it similar in both models. To repeat it. Ml results in 
—J2/J1 = OL > otc when projected onto the spin Hamil- 
tonian (1) and hence in a low-spin GS while M2 resulting 
vcL a < otc and hence in a FM GS. The reason for the re- 
duction of a is mainly a reduction in superexchange due 
to the enhanced Ap^- value. As is clearly seen from Fig. 3, 
whether a < a^ ox a > 0.^ can directly experimentally be 
decided by low-T studies of (j{u)) and such measurements 
are strongly encouraged. In this context we stress once 
more that our main result is not simply the determination 
of a new, improved, parameter set M2 in a literal sense, 
but the demonstration of an unique correlation between 
the magnetic nature of the GS and the T-dependence of 
<j{u>). A refinement of M2 when low-T and H ^ data 
will be available, a generalization of the Hamiltonian (2) 
to include further orbitals or interactions, its application 
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to other cupratcs like Li2ZrCu04 with a different set will 
not change this fundamental interrelation. 

To compare the Ji with those from another 
method, density functional theory (DFT) calculations 
for Li2Cu02 using the structural data of Ref. 5 were 
performed in the local spin density plus Hubbard U 
(LSGA+U) approach. The full-potential local-orbital 
(FPLO) codoi^ was employed. The onsitc Coulomb in- 
tegral for a Cu 3dxy hole was taken to be C/ = 8 eV and 
an onsitc exchange integral J = 1 eV was used. The 
LSDA+[/ yields quite reliable results for total energy 
differences of different spin structures from which the Ji- 
values of Eq. (1) can be extracted. TAB. I shows the Ji- 
and J2-values we found in comparison with similar results 
given in Refs. I9lll4l inch a quantum chemistry study of 
di- and trimers, with those from projecting the models 
Ml and M2 onto Eq. (1). The small Ja-values justify a 
description of Li2Cu02 in terms of a J1-J2 model. There 
is an almost perfect agreement of three independent ap- 
proaches to be close to Uc = 0.25 ( J3 = 0) of a quasi-lD 
situation. Moreover, model M2, the LSDA-|-t/ as well 
as the QC calculations point to a < Ofc, i-e. to a FM 
chain. This holds in our LSDA-hC/, if 7 eV < t/ < 9 eV. 

To conclude, we have shown that the model many- 
body Hamiltonian (2) with parameters M2 is capable of 
consistently describing XAS, EELS, RIXS and o-(w) of 



TABLE I: In-chain exchange integrals Ji and frustration ratio 
a = — J2/J1 obtained from exact diagonalizations (ED) of 
the Hubbard model (2) with sets M2 and Ml as well as from 
independent DFT+Z7 and a quantum chemistry (QC) study. 
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Ji,K 


-146 
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-126 
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J2,K 


33 


30 


49 


52 


22 


J3, K 


-0.5 
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-2 
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0.23 


0.19 


0.47 


0.41 


0.15 



" M2, CU4O10 chain cluster 

*> LSDA-|-i7 (FPLO), (7 = 8 eV, J = 

'^ Ml [1]; CU4O10 chain cluster 

■^ GGK+U (WIEN2fc); U^a = U- J 

= TAB. II of Rcf. 9 



1 eV 



: 8 eV, see Ref. 14 



Li2Cu02 at 300 K. The main issue, however, is that the 
magnetic GS and the spin excitation spectrum strongly 
affect the T-dependence of cr(a;) in the visible range. The 
reason for is the thermal population of excited spin states 
which differ magnetically much from the GS. This qual- 
itatively new effect, irrespectively on details of the rele- 
vant microscopic parameters for a particular compound, 
has been overlooked so far in the cupratc optics liter- 
ature. Our findings allow a quick qualitative magnetic 
classification of cupratcs by optical measurements: if the 
AFM (FM) exchange is dominant, (j{uj) in the ZRS en- 
ergy region increases (decreases) lowering T. In the ID 
FM-J1-AFM-J2 case it even allows a sharp determina- 
tion of the position of a system relative to the critical 



point. Thus, a T-depcndent cr{uj) owing to thermally 
activated spin excitations is a general phenomenon. Ex- 
perimental studies of this T-dependence together with 
a quantitative theoretical analysis of many-body models 
can provide deep insight in the GS and the spin excitation 
spectrum of undopcd insulators. This is relevant also for 
weakly coupled systems with un-sliared Cu04-plaquettes 
such as BisCuOzf^ or Sr2Cu(PO)42i. For LiaCuOa such 
studies can be decisive to settle the intra-chain exchange 
and the reason for the FM in-chain order. li ac > a 
would be true and the FM in-chain order below T/v would 
be due to weak inter-chain exchange, a non-monotonic 
behavior of a{uj,T) above Tjy should be observed. For 
La6Ca8Cu2404i and Y2Ca2Cu50io also showing FM in- 
chain order a similar situation can be anticipated. Near 
a critical point as in Li2ZrCu04i^ also the magnetic field 
dependence becomes an issue. 

We thank the DFG (Drc-269/3-l;D,M), its Emmy- 
Nocther program (R), the ASCR project AVOZ10100520 
(M), the GIF (M,R,E) for support, and M. Knupfer, 
R. Kuzian, J. Richter, and S. Nishimoto for discussions. 

Note added. After submission of the present paper 
a{uj,T) for Sr2Cu03 has been studied by Kim et al. 
(arXiv:0806.2083). There (see Fig. 1) a sfight reduction 
of the intensity of the charge transfer transition (CT) 
near 1.8 eV by raising T has been found, similarly to 
what we have shown in Fig. 3. In our approach this be- 
havior is attributed to optical excitations from thermally 
populated low-lying states which exhibit a reduced (7{lu) 
in the mentioned above CT region. For long chains the 
existence of such states (e.g. triplets) independently on 
the value of Ji > follows from the Bethe-ansatz solu- 
tion of the related ID- AFM spin-1/2 Heisenberg model. 
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